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摘要 : 自由 活塞 斯 特 林 热机 取消 了 传统 斯 特 林 热机 的 曲柄 连 杆 结构 ， 具 有 结构 简单 、 高 效率 、 高 可 靠 性 、 寿 命 长 等 诸多 优点 。 
但 目前 的 研究 热点 多 集中 于 自由 活塞 斯 特 林 发 动机 与 制冷 机 ， 对 自由 活塞 斯 特 林 热泵 的 研究 关注 很 少 。 本 文 设计 了 一 台电 驱动 
自由 活塞 斯 特 林 热 泵 ， 首 次 提出 并 采用 了 热机 部 分 反 向 布置 的 结构 ， 彻 底 避 免 了 高 温 端 换 热 器 对 直线 电机 造成 的 不 利 影响 
给 出 了 实验 结果 。 该 热泵 在 -20~50°C 的 泵 热 温 差 下 获得 了 1.4 的 整 机 COP。 利 用 热 声 理论 ， 本 文 对 该 热泵 进行 了 理论 计算 ， 
与 实验 结果 进行 了 对 比 ， 初 步 考 察 了 热泵 在 不 同 泵 热 温 差 、 不 同 输入 电 功 以 及 不 同 运 行 频 率 的 工作 特性 。 
关键 词 : 斯 特 林 热泵 ; 自由 活塞 ， 反 向 布置 ， 直线 电机 ;， 热 声 
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Theoretical and exper imental research on free piston Stirling heat pump 


LI Ke™”, YU Guo-Yao!, ZHANG Yi-Bing’, LIAO Xu’, DAI Wei ,LUO Er-Cang! 


(1. Key Laboratory of Cryogenics, Chinese Academy of Sciences, Beijing, 100190 China; 
2. University of Chinese Academy of Sciences, Beijing, 100049, China; 
3. Lihan Thermoacoustic Technologies Co. Ltd, Shenzhen, 518005, China) 


Abstract: Free-piston Stirling machine has received much attention due to its high efficiency, fuel flexibility, high 
reliability and long lifetime. However, much research has focused on free piston Stirling generator and free piston 
Stirling cryocooler. Little attention has been paid to the free piston Stirling heat pump. In this paper, we have 
designed and manufactured a 1 kW class free piston Stirling heat pump and tested its performance systematically. A 
reversed arrangement has been put into practice, which could ensure the linear motor working at ambient 
temperature. The COP of 1.4 was achieved with the heating temperature at 50 °C and the ambient temperature at 
-20 °C. Thermoacoustic theory was used for design and calculation and the results agreed well with the experiments. 


Key words: Stirling heat pump; free piston; reversed arrangement; linear motor; thermoacoustic 


0 研究 背景 斯 特 林 热泵 ， 即 用 自由 活塞 斯 特 林 发 动机 驱动 自 
0 活塞 斯 特 林 热 泵 。 而 对 于 采用 直线 电机 了 驱动 的 传统 
自由 活塞 斯 特 林 热 机 是 斯 特 林 热机 的 一 种 具体 结构 自由 活塞 斯 特 林 热 泵 ， 尚 未 见 相 关 报 道 。 针 对 
形式 ， 通 常 分 为 三 种 : 采用 正 斯 特 林 循 环 的 自由 活 这 一 研究 空白 ， 本 文 对 一 台电 驱动 自由 活塞 斯 特 林 
塞 斯 特 林 发 动机 ;采用 逆 斯 特 林 循 环 的 自由 活塞 斯 热泵 进行 了 实验 研究 ， 并 与 热 声 理论 的 计算 结果 进 
特 林 制冷 机 与 热泵 。 自 由 活塞 斯 特 林 热 机 不 仅 具 有 行 了 对 比 。 此 外 ， 为 避免 采用 传统 结构 可 能 引起 直 
传统 斯 特 林 热机 高 理论 效率 ( 卡 诺 效 率 )、 外 燃 性 等 的 线 电机 温度 过 高 的 问题 ， 本 文 首次 提出 并 采用 了 热 
优点 ， 同 时 取消 了 传统 斯 特 林 热机 的 曲柄 连 杆 结构 ， 机 部 分 反 向 布置 的 结构 ， 以 确保 直线 电机 工作 于 室 
解决 了 密封 问题 ， 提 高 了 可 靠 性 及 使 用 寿命 。 因 此 ， 温 环境 。 
自由 活塞 斯 特 林 热 机 一 直 是 国内 外 研究 的 热点 ， 其 、 
中 以 美国 Sunpower 公司 与 Infinia 公司 的 技术 最 为 成 1 日 由 活塞 斯 特 林 发 电机 结构 
熟 且 最 具 代 表 性 。 但 目前 的 研究 多 集中 于 自由 活塞 1(a) 给 出 了 一 种 传统 的 自由 活塞 斯 特 林 热机 
斯 特 林 发 电机 与 制冷 机 ， 对 于 自由 活塞 斯 特 林 热 泵 结构 。 这 种 结构 广泛 应 用 于 发 动机 与 制冷 机 。 对 于 
的 研究 鲜 见 报道 2, 仅 美 国 的 Sunpower 公司 与 日 本 ”发 动机 ，2 为 高 温 端 换 热 器 (又 称 热 头 )，5 为 室温 换 
的 Kawasaki 公司 报道 过 采用 Duplex 结构 的 自由 活塞 热 器 ; 而 对 于 制冷 机 , 2 为 低温 端 换 热 器 (又 称 冷 头 )， 
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5 为 室温 端 换 热 器 。 可 见 ， 无 论 是 作为 发 动机 还 是 制 
冷 机 ， 室 温 端 始终 位 于 直线 电机 侧 ， 从 而 保证 直线 
电机 工作 于 室温 环境 。 但 当 图 1(a) 结 构 直 接应 用 于 热 
泵 时 ， 换 热 器 5 成 为 高 温 端 换 热 器 ， 靠 近 直 线 电 机 ， 
这 会 导致 直线 电机 工作 于 高 温 环 境 ， 不 利于 电机 运 
行 。 针 对 这 一 问题 ， 本 文 提出 了 一 种 新 型 自由 活塞 
斯 特 林 热泵 结构 ， 如 图 1(b) 所 示 ， 将 排出 器 (又 称 自 
由 活塞 ) 反 向 布置 ， 相 应 地 交换 高 温 与 室温 换 热 费 位 
置 ， 使 直线 电机 工作 于 室温 环境 。 图 1 中 还 标示 了 
两 种 热泵 结构 中 的 声 功 流向 与 温度 梯度 ， 具 体 说 明 
见 下 一 节 。 


一 DW Do 


(a) 传 统 结构 


(b) 反 向 布置 结构 
(1. 膨 胀 腔 ，2. 高 温 端 换 热 器 (发 动机 ) 或 低温 端 换 热 器 ( 制 


冷 机 ) 或 室温 端 换 热 器 (热泵 );，3. 回 热 器 4. 排出 器 ;5. 室 温 端 


多 
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21kW 电 驱 动 自由 活塞 斯 特 林 热泵 


Fig.2 1 kW motor driven free piston Stirling heat pump 


Table 1 Parameters of the free piston Stirling heat pump 
表 1 自由 活塞 斯 特 林 热 泵 结构 参数 


部 件 名 称 
室温 端 换 热 器 


二 


高 歇 端 换 热 器 


Ea 


排出 器 


膨胀 腔 


换 热 器 (发 动机 或 制冷 机 ) 或 高 温 端 换 热 器 (热泵 )，6. 压 缩 腔 ; 
7. 动 力 活塞 (或 称 直线 电机 活塞 )，8. 直 线 电 机 ) 

图 1 自由 活塞 斯 特 林 热 机 结构 示意 

Fig.1 Schematics of free piston Stirling machines 

图 2 为 本 文 所 研究 的 自由 活塞 斯 特 林 热 泵 的 实 

物 图 。 该 热泵 设计 频率 为 80 Hz， 工 质 为 氨 气 ， 平均 

压力 5 MPa， 直 线 电 机 最 大 功率 为 1 kW， 各 部 件 具 

体 参 数 见 表 1 。 


现 


压缩 腔 
动力 活塞 
直线 电机 


体 参数 


翅 片 式 紫铜 换 热 器 , 长 35 mm, 200 条 翅 片 ， 
流 道 宽 0.4 mm 
翅 片 式 紫铜 换 热 器 , 长 35 mm, 200 条 翅 片 ， 
流 道 宽 0.5 mm 

长 33 mm， 内 径 33 mm， 外 径 84 mm， 填 
充 180 目 不 锈 钢丝 网 

直径 80 mm, 连接 轴 直 径 25 mm, 动 质量 m 
为 0.52 kg， 板 弹簧 刚度 为 96 kN/m 

45S cc 

28 cc 

直径 60 mm 

动 质量 m 为 1.6 kg， 磁 力 弹 簧 刚度 Kas 为 
59.5 kNm， 背 腔 容积 甩 为 0.9L，BL 数 42 
N/A， 内 阻 R. 为 0.48Q， 电 感到 为 22mH 


2 热 声 理论 


分 析 


近年 来 热 声 理 论 已 经 发 展 成 为 一 种 强大 的 交 变 
流动 分 析 工 具 趾 ， 本 节 将 按照 热 声 理论 , 简要 介绍 自 
由 活塞 斯 特 林 热泵 工作 原理 及 计算 方法 。 


热 声 理论 认为 ， 
部 件 ， 而 并 非 传统 理论 分 析 认 为 的 ， 仅 是 起 到 蓄 热 
与 放 热 作 用 的 热力 海绵 。 热 声 理论 中 的 介 观 热力 学 
理论 外 给 出 了 回 热 器 工作 原理 的 定性 解释 : 众多 气体 


回 热 器 是 实现 功 热 转换 的 核心 


微 团 在 其 平衡 位 置 附近 振动 ， 与 固体 介质 相互 作用 ， 
每 个 气体 微 团 会 经 历 膨胀 、 收 缩 、 吸 热 、 放 热 等 热 
力 过 程 ， 这 些 过 程 形成 了 一 个 个 精巧 的 热力 循环 ， 
实现 了 功 热 转换 ， 再 依靠 气体 微 团 的 接力 作用 ， 最 


终 实现 了 功 热 转换 和 输 运 。 
在 气体 流通 截面 上 对 声 功 强度 进行 积分 ， 
声 


理论 中 一 个 重要 概念 


要 
得 热 


E, =7Rel -| 


即 获 


一 一 声 功 E32， 定义 为 


(1) 


其 中 方 、 
式 ，* 表 示 共 轿 。 声 功 概念 


类 似 传统 热力 学 


| 


示 。 图 1 中 还 给 出 了 发 电机 中 


世 <v 
I 
= 


室温 端 流 向 高 温 端 } 
再 通过 排出 器 传递 至 压缩 腔 ， 在 压缩 腔 中 
声 功 流向 直线 电机 ， 推 动 活塞 做 功 ， 将 


放大 过 程 。 


量 + 一 阶 量 ” 的 形式 ; 

Y= +Re($.e™) 
其 中 ， 了 表示 密度 p、 
理 量 ， 下 标 0 表示 平均 量 前 
动量 部 分 ，w 为 角 频 率 ， 


分 ， 下 标 1 表示 一 
~ 表示 时 域内 波动 量 ， 


程 、 能 量 方程 、 理 想 气体 状态 方程 ， 忽 
阶 以 上 的 小 量 后 可 得 到 如 下 线 
es dp 1 


太 分 别 为 压力 波 与 体积 流 率 的 复数 
分 析 中 的 
PV 功 ， 特 别 地 ， 当 声 功 计算 位 置 位 于 活塞 表面 ， 
得 声 功 等 于 PV 功 。 热 声 理论 还 指出 ， 对 于 发 电机 ， 
回 热 器 中 的 声 功 流向 必须 与 回 热 器 的 温度 梯度 同 
向 ， 才 可 实现 将 热能 转化 为 声 功 ， 如 图 1 中 箭头 所 
\ 体 的 声 功 流向 ， 
在 回 热 嚣 中， 气体 实 现 了 功 热 转换 与 传递 ， 声 功 由 
罕 不 断 增 大 ， 然 后 流入 膨胀 腔 ， 


按照 经 典 热 声 理论 ， 各 个 波动 量 可 写成 “ 平 


温度 T、 压 力 p 与 流速 wx 条 


ROSY 


所 


即 : 


? 大 部 分 
声 功 转化 为 
电能 ， 剩 余 小 部 分 声 功 流向 回 热 器 ， 再 次 进行 声 功 


(2) 
等 物 
阶 波 


示 波 动量 的 复数 幅 值 ， 其 中 包含 着 波动 幅 值 和 相位 。 
将 上 述 简化 形式 的 波动 量 代 入 连续 性 方程 、 动 量 方 
各 二 阶 及 二 
生 热 声 方程 组 (3)~(5) 


© | ©) 
dU d7 ~。 
i 二 cr Ui (4) 
F =2,:P -2,:.0, (5) 
其 中 体积 流 率 U = 站.4，4 为 截面 积 ， 其 具体 参数 
取 值 可 参考 文献 [3]。 
热 声 理论 中 , 另 一 个 重要 的 时 均 量 为 总 能 流 厂 ， 
即 
1 dT 
Hy= FPoCn| Re[Tii]- dA4 (Ak Ak) (6) 
其 中 4 表示 截面 积 ,，k 表示 热 导 率 ， 下 标 f，s 分 别 
表示 流体 与 固体 。 
换 热 器 进出 口 的 总 能 流 之 差 等 于 换 热 器 的 换 热 
量 。 


对 于 排出 器 ， 其 满足 以 下 力学 平衡 方程 ; 
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Dexp Aiisp Peomp (La e A ) I RaispUaisp 


disp 入 (7) 
@ Udisp 二 ne iQOmMaisp laisp 


其 中 天 为 弹簧 刚度 ,Rn 为 机 械 阻 尼 , 下 标 exp、cmp、 
disp 与 rod 分 别 表示 膨胀 腔 、 压 缩 腔 、 排 出 器 与 排 
器 轴 。 


+i 


压 


对 于 直线 电机 可 根据 其 电学 、 力 学 平衡 方程 进 
行 求解 是 。 
U=I(R. +ioL)+ BIA,, (8) 
~ {KV+yp, A 
BLI+i mag b /Po pist ii 
OP (9) 


入 


一 人 List 外 Uist — Po dis = = iQmM ,is pist 
其 中 ，y 为 工 质 气 体 的 比 热 比 ，Zeioaa 为 外 部 电路 负 
载 ， 下 标 pist 表示 电机 活塞 ，masg 表示 磁力 弹 得 。 
具体 的 理论 计算 结果 将 在 第 四 章 中 给 出 。 


3 实验 平台 与 实验 方法 


图 3 为 热泵 实验 平台 示意 图 。 其 中 ， 直 线 电 机 
由 一 台 变 频 电 源 1 驱动 ， 两 台 高 低温 循环 装置 2、3 
用 以 将 传 热 工 质 (水 或 油 ) 维持 在 一 定 温度 ， 实 现 
换 热 器 定 温 换 热 。 泵 热量 O、 输 入 电 功 拨 、 输 入 声 
功 所 、 整 机 COP、 电 机 效率 jmowor、 热 泵 效率 ybp， 
以 及 热泵 相对 卡 诺 效率 ycamo% 可 分 别 按 下 式 计 算 : 


Q=rmc, (Ts, 一 轧 ) (10) 
W.=0.5|i [cos(@, —0,) (11) 
W = 了 pe DO: cos( 9， 人 6) 


(12) 


pist ee i cos(9, -0 + 90") 
20w 2 
COP=Q/W. (13) 
Winotor 三 WW / W. (14) 
m= QO/W, (15) 
NCarnoto% =77 (了 hn 7)/T (16) 


其 中 循环 水 流量 六 由 电磁 流量 计 4 测 得 ， 进 出 口水 
温 TZ,、Tout 由 分 别 由 PT100 铂 电 阻 5、6 测 得 ; 电压 、 
电流 采用 Tektronix 公司 的 电压 传感器 7 与 电流 传 感 
器 8 测量 ， 膨 胀 腔 压力 波动 5, 与 动力 活塞 加 速度 
全 采用 美国 PCB 公司 的 压力 波动 传感器 9 和 加 速度 
计 传 感 器 10 测量 

此 外 ， 在 膨胀 腔 与 压缩 腔 内 分 别 装 有 T 型 热电 
偶 11、12， 用 于 测量 膨胀 腔 与 压缩 腔 内 的 气体 温度 
Tp、Temp。 电 机 背 腔 中 还 安装 有 平均 压力 传感器 13 
用 于 检测 系统 的 平均 压力 po。 

本 热泵 中 采用 了 气 浮 活塞 式 直线 电机 ， 即 依靠 
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气体 轴承 使 活塞 悬浮 于 气 伍 中， 而 气体 轴承 刚度 受 
活塞 位 移 的 影响 。 当 活塞 位 移 过 小 时 ， 和 气体 轴承 不 
足以 支撑 活塞 悬浮 ， 会 造成 电机 活塞 与 气缸 间接 触 
摩擦 。 为 避免 这 一 问题 ， 实 验 中 始终 保证 输入 电 功 
不 低 于 50% 的 额定 功率 。 因 为 未 采用 双 系 统 对 置 结 
构 ， 单 机 的 惯性 力 会 造成 系统 的 强烈 振动 ， 为 此 实 
验 中 将 热泵 固定 于 减 振 平台 ， 如 图 2 所 示 。 

考虑 到 热泵 在 实际 应 用 中 ， 通 常 要 求 使 用 温度 
(高 温 端 温度 ) 固定 ， 而 外 界 环境 温度 会 随 季节 的 
不 同 而 变化 ， 尤 其 是 在 寒冷 地 区 。 因 此 ， 实 验 中 国 
定 高 温 端 温度 为 50*C， 分 别 考察 了 室温 端 温度 、 输 
入 电 功 、 运 


行 运行 频率 对 热泵 性 能 的 影响 。 


(1. 变 频 电 源 ; 2,3. 高 低温 循环 装置 ; 4. 电 磁 流 量 计 ; 5,6. 铂 电阻 ; 
7. 电 压 传感器 ; 8. 电 流传 感 器 ; 9. 压 力 波动 传感器 ; 10. 加 速度 计 ; 
11,12. 热 电 偶 ; 13. 平 均 压力 传感器 ) 


图 3 实验 平台 示意 图 


Fig.3 Schematic of the experimental facility 
PIA E> 士 
4 实验 与 计算 结果 


图 4 给 出 了 在 不 同室 温 端 温度 下 ， 整 机 COP 随 
输入 电 功 的 变化 。 其 中 ， 热 端 温度 维持 在 50*C， 室 
温 端 温度 分 别 为 20、0、-20°C 时 ， 整 机 所 能 达到 的 
最 高 COP 分 别 为 : 1.75、1.5、1.4。 可 以 看 出 ， 实 验 
与 计算 结果 在 变化 趋势 上 基本 一 致 。 对 于 不 同 的 输 
入 电 功 ， 当 室温 端 温 度 为 20°C 时 ， 整 机 COP 随 输 
入 电 功 的 增 大 略 有 减 小 ， 当 室温 端 温 度 为 0 或 -20°C 
时 ， 整 机 COP 基本 不 受 电 功 变化 的 影响 。 对 于 不 同 
的 室温 端 温度 ， 泵 热 温 差 越 大 ， 整 机 COP 越 低 。 
尽管 计算 与 实验 结果 在 变化 趋势 上 保持 一 致 ， 
但 在 具体 数值 上 ， 实 验 中 整 机 COP 仍 比 计算 值 偏 低 
20%~25% 左 右 ， 所 以 有 必要 进一步 对 热泵 与 电机 分 
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别 进行 考察 。 
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图 4 不 同室 温 端 温度 下 整 机 COP 随 输 入 电 功 的 变化 
Fig. 4 COP vs. electric power with different ambient 


temperatures 
图 5 给 出 了 不 同 泵 热 温差 下 电机 效率 随 输入 电 
功 的 变化 。 可 以 看 出 ， 电 机 效率 受 温度 与 电 功 变化 
的 影响 均 不 大 。 但 需要 指出 的 是 ， 实 验 中 仅 考察 了 
电机 在 50% 以 上 额定 功率 下 的 工作 特性 。 当 电机 以 
低 功率 运行 时 ， 如 第 三 节 末尾 所 述 ， 电 机 活塞 与 气 
包间 可 能 产生 接触 摩擦 ， 势 必 造 成 电机 效率 的 降低 。 
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图 5 不 同室 温 端 温度 下 电机 效率 随 输入 电 功 的 变化 


Fig. $ Efficiency of linear motor vs. electric power with different 


ambient temperatures 
此 外 ， 实 验 过 程 中 还 发 现 ， 当 热泵 工作 于 
-20~50°C 泵 热 温 差 时 ， 电 机 机 膏 表 面 出 现 了 明显 的 
结 钉 现象 ， 这 表明 ， 靠 近 电 机 侧 的 室温 端 换 热 器 会 
对 整个 电机 的 工作 温度 产生 强烈 的 影响 。 这 也 从 侧 
面 说 明 ， 对 于 自由 活塞 斯 特 林 热泵 ， 采 用 本 文 所 提 
出 的 反 向 布置 结构 是 十 分 必要 的 。 
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图 6 不同 泵 热 温差 下 热 汞 相对 卡 诺 效率 随 输入 电 功 的 变化 


Fig. 6 Relative Carnot efficiency Vvs. electric power with different 


ambient temperatures 
于 电机 效率 的 计算 与 实验 结果 基本 一 致 ， 所 
以 整 机 COP 的 计算 与 实验 偏差 主要 来 自 热机 部 分 的 
偏差 ， 如 图 6 所 示 ， 这 一 方面 是 由 于 线性 热 声 理论 
的 自身 限制 ， 如 准 一 维 模型 等 ， 另 一 方面 是 由 于 加 
工 工艺 限制 使 热泵 实际 参数 偏离 了 设计 参数 。 不 过 
我 们 仍 能 看 到 实验 与 计算 结果 具有 相同 的 变化 趋 
势 ， 即 相对 卡 诺 效率 随 电 功 变化 较 小 ， 而 受 泵 热 温 
差 的 影响 却 很 大 。 随 着 泵 热 温 差 的 增 大 ， 热 泵 的 相 
对 卡 诺 效率 也 在 升 高 ， 换 而 言 之 ， 自 由 活塞 斯 特 林 
热泵 在 大 泵 热 温 差 下 运行 更 具 优势 。 
图 7 给 出 了 运行 频率 对 于 热泵 性 能 的 影响 ， 其 


也 


工 况 为 室温 端 -200C， 电 功 900 W。 可 以 看 出 ， 实 验 
与 计算 结果 的 变化 趋势 基本 符合 ， 并 表明 ， 对 于 本 
文 所 研究 的 热泵 ， 其 整 机 COP、 电 机 效率 与 相对 卡 
诺 效 率 受 频率 变化 的 影响 很 小 。 
@” 整 机 COP， 实 验 
2.3 e@ 相对 卡 诺 效 率 ， 实 验 100 
4 ”电机 效率 ,实验 。 
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7 整 机 COP、 电 机 效率 、 相 对 卡 诺 效 率 随 频率 的 变化 


Fig. 7 COP, linear motor efficiency and relative Carnot efficiency 


vs. electric power with different ambient temperatures 
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5 结论 


本 文 设计 并 研究 了 一 台 采 用 电 驱 动 的 自由 活塞 
斯 特 林 热泵 ， 首 次 提出 并 采用 了 热机 反 向 布置 结构 ， 
高 温 端 换 热 器 位 于 远离 直线 电机 侧 ， 可 有 效 避 免 高 
温 端 换 热 器 导致 电机 工作 温度 过 高 。 对 该 热泵 在 不 
同 泵 热 温 差 、 输 入 电 功 以 及 运行 频率 下 的 实验 研究 
表明 ， 在 20~50、0~50、-20~50°C 泵 热 温 差 下 整 机 
COP 分 别 约 为 1.75、1.5 与 1.4; 直线 电机 效率 对 于 
泵 热 温 差 、 电 功 及 频率 的 变化 均 不 敏感 ， 热 机 相对 
卡 诺 效率 随 电 功 及 频率 的 变化 很 小 ， 但 受 泵 热 温差 
的 影响 却 较 大 ， 在 20~50、0~50 与 -20~50°C 时 ， 分 
别 为 20%、30% 与 36%。 这 表明 ， 自 由 活塞 斯 特 林 
热泵 可 能 在 高 温 热泵 应 用 领域 更 具 优 势 。 本 文 还 从 
热 声 理 论 的 角度 ， 阅 释 了 该 热泵 的 工作 原理 ， 进 行 
了 理论 计算 ， 并 获得 了 与 实验 结果 相同 的 变化 规律 。 
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